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Abstract. - We study how quantum correlations survive at large scales in spite of their exposition 
to stochastic backgrounds of gravitational waves. We consider Einstein-Podolski-Rosen (EPR) 
correlations built up on the polarizations of photon pairs and evaluate how they are affected by the 
cosmic gravitational wave background (CGWB). We evaluate the quantum decoherence of the EPR 
correlations in terms of a reduction of the violation of the Bell inequality as written by Clauser, 
Home, Shimony and Holt (CHSH). We show that this decoherence remains small and that EPR 
^ ' correlations can in principle survive up to the largest cosmic scales. 
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, It has long ago been suggested by Feynman that gravitation could be at the origin of 

a universal decoherence mechanism preventing macroscopic systems to exhibit quantum 
coherence properties |1]. This idea has been considered since that time through the study 
of different mechanisms [2HZ1, some of them involving Planck-scale physics [5HlO). 
. These important theoretical questions can also be addressed as experimental challenges. 

' In particular the following questions have been discussed in this context [TT] : (i) are 

. quantum interferences limited to microscopic objects having masses or energies smaller than 

some intrinsic limit ? (ii) are quantum correlations limited to experiments performed on 
length scales smaller than some intrinsic limit ? These two questions have led to efforts 
for producing evidence for quantum interferences with larger and larger molecules |12| as 
well as quantum correlations on larger and larger distances |13j . Up to now. the reported 
limits are associated with practical rather than fundamental issues. Efforts are going on 
for mastering these technical limitations and approaching closer and closer the underlying 
fundamental questions in future experiments |1 H I14 ) . 

Here, we will study the effect of the unavoidable interaction with the gravitational waves 
backgrounds which pervade our spacetime environment. It has already been shown that 
such an interaction was responsible for an intrinsic decoherence mechanism in atomic inter- 
ferometry. The mechanism is more and more efficient for larger and larger masses, but it 
depends also on other parameters. With this model of spacetime fluctuations, it is possi- 
ble to associate quantitative estimations to the qualitative Feynman argument, using only 
known gravitational physics, astrophysics and cosmology |15| . 

In the present paper, we give quantitative estimations for the case of Einstein-Podolski- 
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Rosen (EPR) quantum correlations between polarization entangled pairs of photons |16) . 
Precisely, we study how the violation of the Bell inequality pT| , under the form written by 
Clauser, Horne, Shimony and Holt (CHSH) [TH], is affected after the propagation of photons 
over some large distance. We evaluate the effect of the binary confusion background (BCB) 
and of the cosmic gravitational wave background (CGWB) [T^. The former is a classical 
background resulting from the confusion of gravitational waves emitted in the Galaxy and its 
vicinity. The latter, predicted to have been created from quantum fluctuations of the metric 
in primordial cosmology |20| . is essentially characterized by the dimensionless parameter 
rjgw which measures the energy density of gravitational waves with respect to the critical 
density energy. 

The main result of the paper will be that the reduction of the violation of CHSH in- 
equality is bound by the value of fig^ and therefore remains small even after propagation 
at the largest cosmic distances. This means that EPR polarization entanglement is pre- 
dicted to survive at these extreme scales without being washed out by the interaction with 
gravitational waves backgrounds. 

Polarization entangled photon pairs. — The Bell test experiment considered in 
this paper is schematized on the spacetime diagram of fig. [1] A source S (placed for example 
on Earth) sends pairs of polarization entangled photons towards two detectors A and B. The 
geodesic motion of S, A and B is represented by the worldlines which are nearly vertical of 
fig. [T] (their velocities are much smaller than c). The lines inclined at ±45° with respect to 
the space and time axis correspond to the propagation of the photons over a time of flight 
r. 
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Fig. 1: Spacetime diagram associated witii a Bell test experiment with polarization entangled 
photon pairs; the nearly vertical lines represent the motions of the source S and detectors A and B 
while the inclined lines correspond to photons from source to detectors; r is the time of flight from 
S to A or B and T is the integration time needed to estimate correlations. 

The source S emits pairs of polarization entangled photons. We consider that the cor- 
responding two-photon state has a null momentum and a null spin (singlet state). From 
conservation and parity arguments |21| . it is thus deduced that the photon pairs are de- 
scribed by a Bell state 

l^s)-^(|i?^,i?^)-|i^,i^)) , (1) 

where R and L denote right-handed and left-handed circular polarizations while and — > 
refer to the photons propagating from S towards respectively A and B. One of the latter 
observers, say A, measures the polarization of photons with polarization beam splitters and 
photodetectors. He attributes -|- and — values to detections in the two polarization channels, 
when the polarization measurement setup has a local orientation defined by the angle 0a- 
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The origin of this orientation is defined by a local reference. The observers A and B count 

the coincidences over a measurement time T for the 4 possibilities H — h, H — , — h and , 

and they measure probabilities (average numbers for each possibility divided by the number 
of all coincidences) 

sin^e 
cos^ 9 

e = Ob - Oa . (2) 

The two observers A and B then evaluate correlation functions E and S defined respec- 
tively for 2 and 4 orientations of the polarization measurement setups |18) 

^AS. ^ p^^ ^ p _ p^_ _ p_^ ^ _ (-20) ^ 

S^E^^- E^^' + E^'^ + E^'^' . (3) 

Should the conditions of Bell theorem be obeyed [T7] , the following CHSH inequality would 
be deduced [IB] 

\S\ < 2 (4) 

As is well known, this inequality is violated by the polarization entangled photon pairs 
described by the preceding equations (l-|3l). The maximum violation is predicted for example 
with the angles Qa = 0, Ob tt/S, 9^ = 7r/4 and 9^ = Stt/S, which lead to |S'max| = 2^2. 

We want to stress at this point that the effects of gravitation have been ignored in the 
discussion up to now. In particular, the orientation of the polarization measurements setups 
at A and B raises no fundamental question in the absence of spacetime curvature but it has 
to be carefully defined in the context of general relativity. We note that the orientation can 
be controlled by using the expression ([3]), allowing the two remote observers to exploit the 
available quantum correlations to get a relation between the angles 9a and 9b- 

Effect of gravitational fields. — We now consider the same problem in general 
relativity, taking into account the effect of gravitational fields, and in particular gravitational 
waves, on polarization entangled photon pairs. To this aim, we have to study the propagation 
of polarized electromagnetic fields in general relativity. The dominant effect is a rotation 
of the polarization of light along the ray [22]. The use of this effect for the purpose of 
detecting gravitational waves has been studied |23[f21] as well as its possible contribution 
to decoherence [52]. Here, we consider the decoherence of the photon pairs due to the 
interaction with the gravitational waves background. 

In order to treat consistently this problem in general relativity, one has to deal not 
only with the polarization of electromagnetic field but also with the orientation of the local 
reference axis at the two observers. In the eikonal approximation which is sufficient for the 
purpose of this letter, both problems are treated by solving the geodesic equation for the 
motion and parallel transport for the polarization of light and orientation of the polarizers. 
It is only by taking both effects into account that one obtains a properly defined observable, 
associated with a gauge invariant expression. In the following, we use the results of detailed 
calculations presented in |26] . 

A proper observable is defined by comparing the rotation of the photon polarization 
along the ray — > from the source S to the detector B and the rotation of local reference 
orientations at the source and detector. At first order in the gravitational perturbation, this 
observable has the following expression 

a^{t) = - (5i/i23 - 92/ii3)cla . (5) 

^ Jet 
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The integral is taken along the unperturbed path of the photons from the time t to the time 
t + T, with T the time of flight and a the time parameter along this path. The notations 
for the metric are the same as in |27j and the expression ([S]) is written for a propagation 
along the direction x^. The angle and the similarly defined are gauge invariant 
observables. 

Considering first the case of stationary gravitational fields, we see that eqs. (|2l3p are 
changed as follows 



P++ = P— = 



P+- = P- 



2 

E^^ = - cos (2(6 - a)) , a ^ - . (6) 

A comparison of ^ with ^ shows that that stationary gravitational fields lead to a relative 
rotation of the orientations at A and B, taking into account the parallel transport of photon 
polarization as well as that of orientation of the polarizers. The structure of the formula 
can also be understood in analogy with an interferometric signal : the rotation of linear po- 
larizations is indeed equivalent to different dephasings e^'" for the two circular polarization 
states. Using this analogy, we will be able in the next sections to use results already known 
for interferometers and apply them to the case of Bell test experiment. 

Up to now we have supposed that the angle a is time independent or at least that is 
varies slowly enough so that it does not affect the measurements during the integration 
time T. In this case, the already discussed trick can still be used : the two observers at A 
and B can in fact compensate for the slow effect of gravitational fields by adjusting their 
local orientation references. It follows that the maximum violation of Bell's inequality is not 
affected, though the orientation of the polarizers may have to be changed in order to attain 
it. Of course, this can no longer be the case when rapid fluctuations cause a blurring of the 
correlation and lead to a decrease of the maximum violation of CHSH inequality. We focus 
the discussion on this effect in the sequel of this letter. 

Effect of stochastic gravitational waves. — We now discuss the effect of stochastic 
backgrounds of gravitational waves, which produce not only drifts but also fluctuations. 

We flrst remind that two stochastic backgrounds are usually studied. The so-called binary 
confusion background (BCB) is the sum of signals of all unknown binaries in our Galaxy and 
its vicinity. Its stochastic nature is only due to the confusion of a large number of poorly 
known sources. In particular, each binary system produces a classical gravitational wave, 
so that the superposition of all contributions remains a classical stochastic background. In 
contrast, the cosmological gravitational wave background (CGWB) stems from quantum 
fluctuations created during the primordial cosmic era and then amplifled by a huge factor 
through their coupling to the evolution of Universe [20] . 

Gravitational waves may be described in linearized general relativity and the fluctua- 
tions of the metric field described by correlation functions (all notations are as in [27j). The 
backgrounds appear stationary for not too long observations. We also use simple assump- 
tions of unpolarized and isotropic backgrounds, which are valid for the CGWB and first 
approximations for the BCB. All information about the backgrounds is thus contained in 
the one-sided spectral density S'gw [19] 

/•oo 

S'gwH=/ dtcxp{-iut){h{t)h{0)) , (7) 



where h is any one of the amplitudes hi2 or (hn — /i22)/2 (or other ones obtained after 
spatial rotations). The spectrum S'gw can equivalently be described as a spectral density 
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of energy, or a number of graviton per mode Tigw (much larger than unity for all modes of 
interest) , 

^^'■^.wM- (8) 
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For the CGWB, the spectrum Sg^ is often written as 

tS'cr 
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(9) 



where Hq is the present day Hubble rate {Hq = d/a where a is the cosmic size parame- 
ter) while the dimensionless quantity figw represents the energy density of the background 
compared to the present value of the critical energy density for closing the Universe. 
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Fig. 2: Spectra Sgw considered in the discussions : the full (blue) line corresponds to the CGWB 
with f2gw = 10^^'' and the dashed (red) line to the BCB (model given in [29|). 

The parameter flg^ might depend on frequency but it is often considered as constant over 
broad frequency intervals. Though the CGWB has not been detected up to now, constraints 
on the value of may be drawn from a variety of observations |19| . In particular, fig„ has 
to be smaller than a few 10"^'* at w ^ 10~^^Hz (cosmic limit), a few 10~® at w ^ lO^^Hz 
(pulsar limit), and a few 10^® at w ~ lO^Hz (laser interferometry limit) [28]. We have 
drawn as the full (blue) line on Figl2]the spectrum S'gw corresponding to the CGWB with 
rigw = 10"^"*. For the sake of comparison, we have shown as the dashed (red) line on the 
same plot the spectrum corresponding in the BCB, with the model given in |29j . These two 
spectra will be used for some computations in the sequel of this letter. 

Using equations ([5][7|), it is now possible to deduce the correlation function Eg^ now 
obtained after a double averaging over the gravitational environment and over the measure- 
ment time 

i?^„B = -(cos(2(e-c.)), (10) 

where the symbols ( ) and ~ represent respectively a trace over the stochastic gravitational 
waves background and an averaging over the measurement time T. In a linearized treatment, 
the variable a has a gaussian distribution so that the correlation function Eg^ can be 
rewritten 

E^^ = C cos (20 - 2{a)^ , (11) 

C = (cos (2<5a)) = exp {-2Aa^) , 

(5q = a — (a) , Aa^ = {{Sa)'^) . 

The interaction with gravitational waves leads to two different effects : the first one is 
a mean rotation angle (a) which, as already discussed, can be compensated by rotating 
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the polarizers. In contrast, the second effect is a net reduction of the correlation function 
due to the fluctuations 5a which is an exponential function of the variance Aa^ of these 
fluctuations. 

For the CGWB, which corresponds to quantum fluctuations, the variance may then be 
written as the following integral over frequencies 



CGWB 



(12) 



The dimensionless apparatus function A characterizes the response of the experiment to 
gravitational waves at frequency lj. It depends on the geometry of the experiment, and in 
particular on the time of flight t. Its calculation is done in a similar way to that of similar 
functions appearing in the study of decoherence in interferometers [15] or of sensitivity to 
gravitational waves of clock synchronization |30| . At the end of this calculation, it may be 
written 



(13) 



/3^[c.,M] = (l + //)(e-'"^(^-^)-l) . 

The amplitudes measure the sensitivity to the gravitational mode at frequency lu and 
wavevector with /i = ck^ juj the direction of (^ is defined for the photon — >, and replaced 
by — /i for the photon ). 

After straightforward integrations, the function A is finally written as 



A{u\ 



5 

16 
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24 + 6cos(2wr) 



15sin(2a;r)\ 



(14) 



The result is drawn on Fig. [3] as a function of cjt. It tends to the constant 5/2 at the limit 
of high frequencies and behaves as (a;r)^/21 at the limit of low frequencies. For the purpose 




Fig. 3: Apparatus function Ar^ drawn versus (ixjt). 



of comparing the two effects, we give also the variance for the BCB, which corresponds to 
classical fluctuations and leads to a slightly different expression 



Aa 



BCB 



(15) 



J^[a;] = 1 — I sine 



sinc(.x) 



smx 
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The two variances ([T2|) and (|15|) are integrals of the gravitational waves background over 
frequencies to larger than the inverse of one of the typical time involved in the correlation 
measurement. For the CGWB, this time is just the time of flight r of the photons from the 
source to the detectors. For the BCB, this time is the longer one of the time of flight r and 
the measurement time T . The effect associated with r may be considered as intrinsic, as 
it cannot be reduced for a given geometry. In contrast, the effect associated with T can in 
principle be limited by having this measurement time shorter than r. 

Discussion and conclusions. — We now come to the discussion of the decrease of EPR 
quantum correlations, i.e. the reduction of the maximal violation of the CSHS inequality 
(|3]), due to the effect of the CGWB. For completeness, we also present results for the BCB. 

From eq. ljlip . we deduce the maximal value S'max of the parameter S*, after optimal 
angles have been chosen, 

S'max = C 2V2 = 2V2 cxp (-2Aa2) . (16) 

Following the analogy between the EPR correlation function and interferences, it turns out 
that the reduction of the maximal violation of Bell-CSHS inequality is analogous to the 
decoherence effect already studied for interferometers [15]. The main difference is that the 
reduction is determined by the variance of an angle a whereas it was determined by the 
variance of a phase for interferometers. 

We show on Fig|4] the results of the evaluation of the quantity 2Aa^ which appears in 
the exponential factor in eq. ([T6| . The full (blue) line corresponds to the effect ([T2|) of the 
CGWB (flgw = 10^1^), and the dashed (red) hue to the effect ^ of the BCB ((model given 
in [29])). All curves are drawn as functions of the time of flight t varying from 10~'^s - typical 
time of flight from a station on Earth to the ISS - to 10*'s - time of flight corresponding 
to the scale of the solar system (sa 2000 astronomical units). For the case of the BCB, the 
various curves correspond to different values of the measurement time T varying also from 
10~'^s to lO^s. As expected all variances are increased when the time of exposition (r or T) 
to the gravitational wave backgrounds are increased. The effect of the BCB is larger than 
that of the CGWB for not too large times of exposition, which just reflects the shapes of 
the corresponding spectra which were plotted on Fig|21 
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Fig. 4: Variances 2Aq^ appearing in the exponential factor in eQ. (|16p plotted as functions of the 
time of flight r : the full (blue) line corresponds to the effect of the CGWB (figw = 10""), 
and the dashed (red) line to the effect (|15p of the BCB ((model given in [29)11 ; in the latter case, 
the various curves correspond to different values of the measurement time T — 10"'^, 1, 10^ and 
lO'^s. 

We see on Fig|4]that the variance Aa^ remains extremely small for any scale accessible 
to experiments in the solar system. This entails that the quantum correlations coded in 
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polarization entangled pairs of photons should survive the exposition to gravitational wave 
backgrounds in any presently foreseeable experiment. In fact, this conclusion can be made 
even more general through an argument which sheds light on an unexpected deep connection 
between the discussion of large scale EPR correlations and cosmology. 

In order to write down this argument, we may rewrite the reduction factor C appearing 
in eas. (jllll6p as follows, under the assumption of a constant value for Jlgw; 

When the form ([H)) of A is used, the integration leads to the closed expression 

T=\{H,rY. (18) 

The resulting r(r) is the straight line appearing on the log-log plot of FiglSl which we have 
drawn for r varying up to the Hubble time t ~ lO^^s. 
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Fig. 5: Dimensionless factor F as a function of the time of flight r on log-log scales; the full (blue) 
line corresponds to the CGWB, while the dotted and dashed (red) lines corresponds to the BOB 
with T < r and T = 1000s respectively. 

FigO shows that the number F is extremely small at all scales small with respect to the 
Hubble scale. In fact, F increases as the square of the time scale of the experiment and it 
would approach unity only at the largest cosmological scales. When evaluating the reduction 
factor C in ([TT]), the small number F has first to be multiplied by a second small number 
r2gw and then exponentiated. As a consequence, the reduction has to remain negligible in 
any experiment. A different conclusion could only be obtained at the largest cosmological 
scales, if the energy of gravitational waves would be a dominant constituent of the content 
of the Universe. 

As a conclusion, we have shown that the quantum correlations coded in polarization 
should survive the exposition to gravitational wave backgrounds. In particular, initial quan- 
tum correlation arising from primordial cosmic processes are still present in current day 
Universe even if their detection is a challenge. Note that these conclusions have been ob- 
tained for gravitational wave backgrounds as they are predicted by standard calculations, 
and they would of course be affected if larger fluctuations were produced by physical pro- 
cesses beyond the standard model. Note also that this conclusion relies on a fascinating and 
unexpected connection between the reduction of EPR quantum correlations at large scales 
and the fundamental cosmic parameter f2gw. 
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